
*dd 2- 3 AUG 2(M 



:u5!ttim'£u (Si-sjiii^'f.miiuiiauisaj Liun^iPTr gngn: tern gnmnm'i tn ir 



Mil 












v • ^ 






eh 







UNITED STATES DEPARTMENT OF COMMERCE 
United States Patent and Trademark Office 

July 22, 2004 

THIS IS TO CERTIFY THAT ANNEXED HERETO IS A TRUE COPY FROM 
THE RECORDS OF THE UNITED STATES PATENT AND TRADEMARK 
OFFICE OF THOSE PAPERS OF THE BELOW IDENTIFIED PATENT 
APPLICATION THAT MET THE REQUIREMENTS TO BE GRANTED A 
FILING DATE UNDER 35 USC 111. 



APPLICATION NUMBER: 60/514,928 
FILING DATE: October 28, 2003 



1 PRIORITY DOCUMENT 

i SUBMITTED OR TRANSMITTED IN 
! COMPLIANCE WITH 

! RULE 17.1(A) OR (B) 




By Authority of the 

COMMISSIONER OF PATENTS AND TRADEMARKS 




e: 

Certifying Officer 



mm 



0 



Please type a plus sign (♦) Inside this box 

PTO/SB/16 (02-01) 
Approved for use through 10/31/2002. OMB 0851-0032 
Patent end Trademark Office; U.S. DEPARTMENT OF COMMERCE 
Under the Paperwork Reduction Act of 1 995. no persons are required to respond to a collection of information unless U displays a valid OMB control number. 

PROVISIONAL APPLICATION FOR PATENT COVER SHEET 
This Is a request for filing a PROVISIONAL APPLICATION FOR PATENT under 37 CFR 1.S3 (c). 



| Express Mail Label No. EL998012521US 



£00 
.OJ 
WO) 

SO 



INVENTORY) 



Given Name (first and middle pfanyjj 



Family Name or Surname 



Residence 
(City and either State or Foreign Country) 



Guy A. 



LAVI 



NORDIA, ISRAEL 



□ Additional inventors are being named on the _ 



. separately numbered sheets attached hereto 



TITLE OF THE INVENTION (280 characters max) 



MAPPING THE CORONARY ARTERIES ON A SPHERE 



Direct aO correspondence to: 
□ Customer Number 
OR 



CORRESPONDENCE ADDRESS 

i — - 



Type Customer Number here 



Place Customer Number 
Bar Code Label here 



ft* Firm or 

Individual Name 



THOMAS E KOCOVSKY, JR 



Address 



FAY. SHARPE, FAGAN. MINNICH & McKEE. LLP 



Address 



1100 SUPERIOR AVENUE. SEVENTH FLOOR 



City 



CLEVELAND 



State 



OH 



ZIP 



44114 



Country 



US 



Telephone 



216/861-5582 



Fax 



216/241-1666 



ENCLOSED APPLICATION PARTS (check all that apply) 



I 12 I ~ ricDtsl. Number ! T 
I 6 I □ Other (specify) l I 



| Specification Number of Pages 
| Drawing(s) Number of Sheets 



O Application Date Sheet See 37 CFR 1.76 



METHOD OF PAYMENT OF FILING FEES FOR THIS PROVISIONAL APPLICATION FOR PATENT (check one) 



n Applicant claims small entity status. See 37 CFR 1.27. 
A check or money order is enclosed to cover the firing fees 

[>3 The Commissioner is hereby authorized to charge filing r- 
fees or credit any overpayment to Deposit Account Number L 
n Payment by credit card. Form PTO-2038 is attached. 



FILING FEE 
AMOUNT ($) 



OrVOSOfl 



JLfiQl 



The invention was made by an agency of the United States Government or under a contract with an agency of 
the United States Government 



□ No. 



□ Yes, the name of the U.S. Government agency and 




TYPED oTPRlNTED NAME 
TELEPHONE 216 861-5582 



ocket Number: j PHUS030268USp] 



USE ONLY FOR FILING A PROVISIONAL APPLICATION FOR PATENT 



This ccflectten of Information Is required by 37 CFR 1.51. The Ertfonnafinn ts used by the pupG c to fila (and by fits PTO to process) a provisions] nppncaflfln. 
ConftSenlfe&tyte governed by 35 US.C. 122 and 37 CFR 1.14. Trte ccCecttoo Is es&nated to tate 8 hours to compete, tndudftg gathartng, preparing, and sufcmitflng 
tho complete provisionaJ application to the PTO. Tlma wtB vary depending upon tta individual case. Any comments on the amount of 8ma you require to complete this 
farm ond/or suggestions for reducing this burden, should be sent to the Chief tafa nn e flo n Officer. U.S. Patent end Trademark Office, U.S. Department of Commerce, 
Washington, 0.C, 20231. DO NOT SEND FEES OR COMPLETED FORMS TO THIS ADDRESS. SEND TO: Box Provisional Appflcaficn. Assistant Ottnmlsstoner for 
Patents. Washington, D.C 20231. 



Page 1 of 1 



Attorney's Docket No. PHOS030268OSP (#899946) 



APPLICATION DATA SHEET 

Mail Stop: Provisional Patent Application 
Commissioner for Patents 
P.O. Box 1450 

Alexandria, VA 22313-1450 
Inventor Information 



Inventor One Given Name : : 
Family Name : : 
City: : 
Country: : 

Correspondence Information 

Name Line. One : : 
Name Line Two : : 
Name Line Three : : 
Address Line One:: 
City: : 

State or Province:: 
Postal or Zip Code:: 
Telephone One:: 
Fax: : 

Electronic Mail:: 

Representative Information 

Registration Number: : 

Application Information ' 

Title Line One : : 
Total Drawing Sheets : : 
Application Type:: 
Docket Number: : 



Guy A, 
LAVI 
Nordia 
ISRAEL 



Thomas E. Kocovsky, Esq. 

FAY, SHARPE, FAGAN, 

MINNICH & McKEE, LLP 

1100 Superior Avenue, Seventh Floor 

Cleveland 

OH 

44114-2518 

(216) 861-5582 

(216) 241-1666 
TKocovsky@faysharpe. com 



28,383 



MAPPING THE CORONARY ARTERIES ON A SPHERE 
6 

Provisional 

PHUS030268USP (899946) 



"Express Mail" Mailing Label Number EL998012521US 

Date of Deposit: OCTOBER 28. 2003 

I hereby certify that this paper or fee is being deposited 
with the United States Postal Service "Express Mail 
Post Office to Addressee" service under 37 C.F.R. 1.10 on the 
date indicated above and is addressed to: Provisional Patent Application, 
Commissioner For Patents. P.O., Box 1450, Alexandria, VA 22313-1450. 



By: Hilary M. McNULTY U 



T>! \nr.VPF.K\PKR2B22 . DAT 



PHUS030268USP 



MAPPING THE CORONARY ARTERIES ON A SPHERE 

Background of the Invention 

The present invention relates to the anatomical imaging arts. It 
finds application in conjunction with the display of the entire coronary tree in the 
heart context and will be described with particular reference thereto. Of course, the 

5 technique is also applicable to mapping other anatomical regions, such as the 
cerebral cortex. Although described with particular reference to CT tomography, it 
will further be appreciated that the invention is equally applicable to other 
diagnostic imaging techniques which generate two and three-dimensional digital 
diagnostic images including coronary arteries for analysis. 

10 A variety of imaging modalities can be used to determine structural, 

functional and perfusion characteristics of the heart and coronary arteries. 
Substantial research and development have been concentrated on new technologies 
and methods that can aid the diagnosis and analysis of the coronary arteries, and on 
planning appropriate treatment. Much of this effort has been focused on the 

15 improvement of the visualization of the coronary arteries and the development of 
computer-assisted analysis of the resulting images. 

Imaging and analysis of the coronary arteries serve several 
purposes: visualization of the coronary tree, detection and quantification of 
stenosis, quantification of the coronary tree reserve, and analysis of the vessel 

20 walls. Traditionally, x-ray angiography has been used to detect and quantify 
stenosis in the coronary arteries. The development of multi-slice CT scanner 
technology having increasingly improved resolution particularly in the slice 
direction made the CT imaging of vascular systems attractive for clinical 
applications. Multislice detectors allow several slices to be acquired and processed 

25 simultaneously enabling more accurate quantification of stenosis. 

Many visualization techniques and quantification strategies are 
currently in use. One approach is the construction of multi-planar reformatted 
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images (MPR) and/or slab maximum intensity projections (slab-MIP) in optimally 
chosen planes, so that the major segments of the coronary tree are visualized. In 
addition, the CT's ability of acquiring 3D volumetric images can be used to 
acquire volume rendered images for visualizing an isolated coronary arteries tree. 

5 However, the current methods for coronary arteries inspection, e.g. curved MPR 
visualization of a single vessel, the planar slab MIP visualization and the volume- 
rendered visualization of an isolated tree either lack the completeness or the 
context of the coronary tree. 

There is a need for a visualization technique that will allow the 

10 volumetric display of the entire coronary arteries tree in context, i.e., location, 
connectivity and surroundings. The present invention provides a new imaging 
apparatus and method which overcomes the above-referenced problems and others. 

i 

Summary of the Invention 

IS In accordance with one aspect of the present invention a diagnostic 

imaging system for displaying a vessel tree is disclosed. A means defines a base 
surface. A means grids the base surface to define pixels. A means projects along a 
normal of each pixel. A means assigns each pixel a grayscale value based on 
grayscale value of voxels intersected by a corresponding normal. 

20 In accordance with another aspect of the present invention a method 

of displaying the coronary arteries tree is disclosed. A base surface is defined. The 
base surface is gridded to define pixels. A normal of each pixel is projected along. 
Each pixel is assigned a grayscale value based on grayscale value of voxels an 
associated normal intersected. A true surface is determined. 

25 One advantage of the present invention resides in displaying the 

entire coronary arteries tree in its context. 

Another advantage resides in using a closed non-planar surface as a 
base for displaying the data of interest 

Another advantage resides in minimization of distortions in the 

30 visualization modes. * 

Another advantage resides in enabling the user to explore the entire 
structure of the cardiovascular system at once. 
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Another advantage resides in permitting the user to adjust the 
thickness of the slab MIP to include more or less data in the image presented. 

Still further advantages and benefits of the present invention will 
become apparent to those of ordinary skill in the art upon reading and 
5 understanding the following detailed description of the preferred embodiments. 

Brief Description of the Drawings 

The invention may take form in various components and 
arrangements of components, and in various steps and arrangements of steps. The 
10 drawings are only for purposes of illustrating the preferred embodiments and are 
not be construed as limiting the invention. 

FIGURE 1 is a diagrammatic illustration of a diagnostic imaging 
system in accordance with the present invention; 

FIGURE 2 is a diagrammatic illustration of a portion of a diagnostic 
15 j imaging system in accordance with the present invention; 

FIGURE 3 is an illustration of a base sphere gridded in accordance 
with the present invention and arteries in the vicinity of the sphere; 

FIGURE 4 is an illustration of a coronary arteries tree in a true 
amorphous surface form; 
20 FIGURE 5 is an illustration of a coronary arteries tree in a globe 

surface form; and 

FIGURE 6 is a 2D map of a coronary arteries tree. 

Detailed Description of the Preferred Embodiments 

25 With reference to FIGURE 1, a diagnostic imaging apparatus 10 

generates electronic diagnostic image representations. In the preferred 
embodiment, the diagnostic imaging apparatus 10 includes a CT scanner which 
generates a three-dimensional volumetric image representation which is made up 
of voxels whose position is defined along orthogonal axes. Other diagnostic 

30 imaging systems, such as magnetic resonance imaging, PET imaging, SPECT 
imaging, and other diagnostic techniques which generate analogous three- 
dimensional diagnostic images are also contemplated. 



In the preferred embodiment, the CT scanner includes a stationary 
gantry 12, in which a rotating gantry 14 is mounted. The rotating gantry carries an 
x-ray tube 16 and a two-dimensional array of detectors 18, which are diametrically 
opposed to each other across a scan circle or an examination region 20. A subject 
5 support 22 supports a region of interest of the subject in the examination region 20. 
Prior to imaging, the subject is preferably injected with a known contrast agent that 
produces the voxels of the known intensity in the vessels, preferably, at the 
intensity higher than that of tissue. A longitudinal drive motor 24 moves the 
subject support 22 longitudinally through the examination region 20. Preferably, 
10 the subject support 22 longitudinally reciprocates the subject as the rotating gantry 
14 rotates continuously for spiral scanning. The imaging apparatus generates a 
volumetric image representation of a transverse volume of the subject, which 
includes the subject's heart, or other region of interest Alternately, the subject 
support 22 can be stepped and data can be collected along a series of parallel, 
15 transverse slices. Although the detector array 18 is illustrated as rotating with the 
rotating gantry 14, the detector array 18 may optionally be mounted as a 
continuous ring on the stationary gantry 12. 

A three-dimensional reconstruction processor 26 reconstructs the 
output signals from the detector array 18 in accordance with the angular position of 
20 each detector element, the angular position of the x-ray tube 16, and the 
longitudinal position of the subject support 22 at the time of sampling into one or 
more three-dimensional image representations. The reconstructed image 
representation is stored in a volumetric image memory 28. 

With continuing reference to FIGURE 1 and further reference to 
25 FIGURES 2 and 3, a base surface means or processor 30 approximates a base 
surface 32, which is preferably a sphere. Of course, other surfaces such as 
ellipsoids, anatomical heart approximations, and the like, are also contemplated. 
i The approximation is based on the assumption that the coronary arteries lie on a 
substantially smooth closed surface. More specifically, a volume selecting means 
30 34 selects a volume of data from the volume memory 28 which corresponds to the 
region of interest, e.g. heart, brain, or other organ. A centerlines determining 
means 36 finds the centerlines 38 of the vessels in the selected volume by one of 
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. known techniques, A best fitting surface means or process 40 draws a best fitted 
surface, preferably the sphere or ellipsoid, to the determined centerlines. Some of 
vessels will lie above the surface of the base surface 32, some of the vessels will 
lie underneath it, and some will have a portion above and a portion below the 
5 surface. The sphere or ellipsoid is rotated such that the axis of rotation is 
substantially parallel to a long axis of the left ventricle- In performing the best 
fitting process 40, a centerlines coordinates converting means 42 converts the 
centerlines coordinates to spherical coordinates according to equations 1. 

10 p = atan[Z/V(X 2 + Y 2 )] Equations 1 

X-a tan [Y / X] 
A = [V(X 2 + Y 2 )/cosp]-R 

where <p is the latitude; X is the longitude; h is the distance from the sphere; X, Y, 
15 Z are the Cartesian coordinates of a centerline point; and R is the radius of the 
sphere. 

With continuing reference to FIGURES 1, 2 and 3, a gridding 
means 50 spreads a grid over the base surface 32 to grid the sphere, or other 
determined surface, into pixels 52. The base surface 32 is stretched or shrunk 

20 along the sphere normals to fit the true form of the vessels wherever a centerline is 
located such that the stretched surface lies on the centerlines. With further 
reference to FIGURE 4, the stretching and shrinking of the sphere to fit the 
centerlines is performed by a true surface determining means 54 which builds a 
true surface 56 representing a "True Form'* undistorted mode visualization, in 

25 which the coronary arteries tree in its context is revealed over a true form surface 
running through its vessels. More specifically, a surface adjusting means 58 
receives the centerlines information from the centerline means 36 and moves the 
location of the spherical base surface 32 along the normals, perpendicularly to the 
surface, to the point where the normals intersected associated centerlines of the 

30 vessels. The remainder of the spherical base surface is mapped in and out to define 
the smoothly varying true form surface. The resulting true surface is the arbitrary, 
amorphous surface, as seen in FIGURE 4. 
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A thickness determining means 60 determines a wall thickness of 
the true surface 56. Preferably, a user supplies predetermined thickness value to 
the thickness determining means 60, e.g. the thickness of the hollow true form 
volume. Alternatively, the thickness is determined by the thickness determining 

5 means 60 based on specified criteria, such as viessel diameter, or other 
characteristics of the subject. A normal projecting or tracking means 62 projects a 
normal from each pixel 52 orthogonally with respect to the true surface 56 in both 
directions, e.g. outward and inward from the true surface 56, to the limits of the 
wall thickness determined by the thickness determining means 60. A grayscale 

10 assigning means 70 examines the intensity of each voxel along the normal and 
assigns a grayscale or intensity value to the corresponding pixel of the true surface 
56 to be displayed in the resultant image. A screening means 72 screens the 
grayscale values of voxels, intersected by each normal, based on predetennined 
criteria to select a grayscale intensity value which will be displayed for the 

15 corresponding pixel. In the preferred embodiment, the screening means 72 selects 
a maximum intensity value along each normal and stores it in a maximum 
intensity image pixel memory 74. The cycle repeats until all the normals are 
screened and a grayscale value of each pixel is determined. 

With continuing reference to FIGURE 1 and further reference to 

20 FIGURE 5, a globe means 80 presents a "Globe" mode of visualization, in which 
the coronary arteries tree is shown on a sphere which is best fitted to the 
determined amorphous true surface 56. More specifically, a draping means 82 
drapes or projects the grayscale values, which have been assigned by the grayscale 
assigning means 70 to the true form surface onto the spherical base surface by 

25 projecting the voxels along the normal to the sphere corresponding to each 
spherical or globe surface pixel. The resulting surface is a globe like structure 84 
and may be rotatably examined as a globe. 

To display the 'True Form" or "Globe" on a conventional monitor 
90 having a two-dimensional display, a video processor 92 projects a normal from 

30 each pixel of the display. Each display pixel is assigned the grayscale value of the 
true surface or the globe surface that it intersects. It will be noted that surface 
pixels closest to the viewing screen are displayed substantially undistorted, while 
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surface pixels toward the poles and the periphery become progressively more 
compressed to give the viewer the impression of viewing a 3D object on the 2D 
viewing screen. Optionally, the 3D effect can be enhanced by assigning an 
illumination direction and adding shading. 

5 With continuing reference to FIGURE 1 and further reference to 

FIGURE 6, a 2D means 100 maps the globe 84 or a selected portion of it on a flat 
surface analogous to the way that the round surface of the earth is projected onto 
flat maps. More specifically, a 2D gridding means 102 spreads the original grid of 
the base surface 32 on a 2D x-y surface. Typically, equatorial surface pixels are 

10 substantially undistorted, while more polar surface pixels are stretched to maintain 
spatial relationships. A matching means 104 matches each coordinate of the base 
surface 32 to corresponding coordinates on the x-y viewing plane. A 2D grayscale 
processor or means 106 assigns each pixel on the x-y viewing surface a grayscale 
value of the corresponding pixel from the globe surface 84. 

15 A mode selection means 108 receives an input from a user and 

selects a visualization mode for processing. The selected visualized data is 
processed by the video processor 92 and displayed then on a monitor 90. A user 
input device 110, such as a keyboard, enables the user to control the video 
processor 92 to display selected projections, rotate the globe or true surface, and 

20 input necessary information as was discussed above. 

In one alternate embodiment, the true form surface and its wall 
thickness are defined, but then data is projected directly onto the spherical (or 
elliptical) base surface. In another alternate embodiment, the two dimensional 
image is mapped directly from the true form image rather than the globe image, 

25 which was derived from the true form image. 

The invention has been described with reference to the preferred 
embodiments. Modifications and alterations will occur to others upon a reading 
and understanding of the preceding detailed description. It is intended that the 
invention be construed as including all such modifications and alterations insofar 

30 as they come within the scope of the appended claims or the equivalents thereof. 
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CLAIMS: 

Having thus described the preferred embodiments, the invention is 
now claimed to be: 

1. A diagnostic imaging system for displaying a vessel tree 

comprising: 

a means (30) for defining a base surface (32); 
a means (50) for gridding the base surface to define pixels (52); 
a means (62) for projecting along a normal of each pixel; 
a means (70) for assigning each pixel a grayscale value based on 
grayscale value of voxels intersected by a corresponding normal. 

2. The system as set forth in claim 1, wherein the base surface 
defining means (30) includes a means (36) for determining vessels centerlines (38) 
and further including: 

i 

a means (54) for mapping the base surface (32) to the centerlines 
(38) to define a true form surface (56). 

3. The system as set forth in claim 2, further including: 

a means (60) for defining a wall thickness to the true form surface 

(56). 

4. The system as set forth in claim 3, wherein the grayscale 
assigning means (70) assigns each pixel (52) a maximum of grayscale values of 
voxels within the defined wall thickness intersected by the corresponding normal. 

5. The system as set forth in claim 2, further including: 

a means (80) for determining a globe surface (84) including a means 
(82) for mapping the assigned grayscale values into a spherical surface. 0 



6. The system as set forth in claim 5, further including: 

a means (100) for projecting the globe surface (84) into a two 
dimensional surface. 

7. The system as set forth in claim 6, wherein the projecting 
means (100) includes: 

a matching means (104) which matches coordinates of the spherical 
surface to coordinates of the two dimensional surface; and 

2D grayscale processor (106) which assigns each pixel on the two 
dimensional surface a grayscale value assigned to at least one corresponding pixel 
on the globe surface (84). 

8. The system as set forth in claim 7, further including: 

a means (108) for selecting at least one of the true form surface, the 
globe surface and the two-dimensional surface for displaying on a monitor (90). 

9. The system as set forth in claim 1, wherein the base surface 
(32) is a sphere or ellipsoid. 

10. A diagnostic imaging apparatus (10) comprising: 

a scanner which examines a region of a subject including coronary 
arteries and acquires three-dimensional data; 

a reconstruction processor for reconstructing the three-dimensional 
image data into a volumetric three-dimensional image representation; 

the diagnostic imaging system of claim 1 for converting a portion of 
the three dimensional image representation into a coronary arteries tree display; 
and 

a display (114) connected to the diagnostic imaging system of claim 
1 for displaying the coronary arteries tree in a context of the region of interest. 
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11. A method of displaying the coronary arteries tree 

comprising: 

defining a base surface; 

gridding the base surface to define pixels; 

projecting along a normal of each pixel; 

assigning each pixel a grayscale value based on grayscale value of 
voxels an associated normal intersected; and 
determining a true surface; 

12. The method as set forth in claim 11, wherein the step of 
defining the base surface includes: 

obtaining a substantially spherical volume data; 
determining locations of centerlines of vessels in the volume data 
based on predetermined grayscale value; and 

generating a best fitted surface through the centerlines. 

13. The method as set forth in claim 12, further including: 
translating the base surface along the normals to overlie points, in 

which the normals intersected associated centerlines; and 

defining a spherical thickness which extends in both directions of a 
boundary of the translated surface. 

14. The method as set forth in claim 13, further including: 
injecting a subject with a known contrast agent which produces the 

highest intensity value inside the vessels; 

in the step of projecting, searching for points with the highest 
intensity the associated normal intersected in the determined thickness; and 

assigning each pixel a maximum intensity value chosen from a 
plurality of grayscale values of voxels the associated normal intersected in the 
determined thickness. 



15. The method as set forth in claim 14, further including: 
mapping the determined maximum intensity values into the 
translated surface. 



16. The method as set forth in claim 14, further including: 
draping the assigned maximum intensity values into the base 

surface to create a globe image; and 

displaying the coronary arteries tree in the globe image which is 
rotatably visualized. 

17. The method as set forth in claim 11, wherein the base 
surface is a sphere. 

18. The method as set forth in claim 11, wherein the base 
surface is an ellipsoid. 

19. A scanner for diagnostic imaging including: 

a stationary gantry which defines a subject receiving aperture; 

a source of an x-ray radiation rotatably mounted on the gantry, 
which source transmits x-ray radiation through a subject disposed in a 
subject receiving aperture; 

a two-dimensional radiation detector for detecting radiation 
transmitted by the source after passage of the radiation through the subject in 
the subject receiving aperture; 

a reconstruction processor which reconstructs x-ray radiation 
received by the two-dimensional radiation detector into a volumetric image 
representation; 

an image processor for performing the method of claim 11; and 
a display for displaying the coronary arteries tree in a context of 
the region of interest. 
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MAPPING THE CORONARY ARTERIES ON A SPHERE 

Abstract 

A coronary arteries tree is approximated by a base sphere (32) 
which is best fitted to vessels centerlines (38). The base surface (32) is gridded to 
define pixels (52). The base sphere (32) is mapped to fit the centerlines (38) such 
5 that a true form surface (56) is determined. A wall thickness to the true form 
surface (56) is defined, preferably, by a user. A normal of each pixel (52) is 
searched for grayscale values of voxels. Each pixel (52) is assigned a maximum of 
grayscale values of voxels within the defined wall thickness intersected by the 
corresponding normal. The resulting true form surface is undistorted mode of 
10 visualization revealing the arteries tree in its context running on the true surface 
drawn through the vessels. Mapping the assigned grayscale values onto the base 
sphere (32) visualizes arteries tree on a globe surface (84) which might be rotatably 
inspected as a globe. Mapping the assigned grayscale values into a flat surface 
visualizes arteries tree on a two-dimensional map. 
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